KL-TH-95/16 



bep-th/9510046 



BRST-INVARIANT 
APPROACH TO QUANTUM MECHANICAL TUNNELING 



JIAN-GE ZHOU, F. ZIMMERSCHIED, J.-Q.LIANG*and H.J.W. MULLER-KIRSTEN 

' Department of Physics, University of Kaiserslautern, P.O.Box 3049 

^\ \ 67653 Kaiserslautern, Germany 

ON ■ E-mail: zimmers®. physik.uni-kl.de 



> 



in 

ON 



X 



ABSTRACT 



O 

ON ' A new approach with BRST invariance is suggested to cure the degeneracy 

problem of ill defined path integrals in the path-integral calculation of quantum 
mechanical tunneling effects in which the problem arises due to the occurrence 
of zero modes. The Faddeev-Popov procedure is avoided and the integral over 
the zero mode is transformed in a systematic way into a well defined integral 
CD ■ over instanton positions. No special procedure has to be adopted as in the 

Faddeev-Popov method in calculating the Jacobian of the transformation. The 
quantum mechanical tunneling for the Sine-Gordon potential is used as a test 
of the method and the width of the lowest energy band is obtained in exact 
agreement with that of WKB calculations. 



1. Introduction 



Quantum tunneling has attracted considerable interest because of its wide appli- 
cation in areas ranging from condensed matter to high energy physics. The instanton 
method is a powerful tool for the calculation of tunneling effects. Recent interests in 
tunneling effects were initiated by the work of RingwaldB, who argued that although 
the cross section of the standard electroweak theory is proportional to an exponen- 
tially small WKB suppression factor, it is nevertheless rapidly growing with energy 
due to multiple production of Higgs and vector bosons. However, it has been observed 
that instantons interpolate between neighbouring vacua and satisfy vacuum boundary 
conditions and therefore may not be adequate for a description of tunneling at high 
energyn. Motivated by the instanton-induced baryon-number violating processes the 
instanton method of quantum tunneling has recently been extended to-,tunneling at 
finite, nonzero energy by means of nonvacuum or periodic instantonsB which have 
become the subject of extensive investigation under the name of sphaleronsufl. 

It is well known that due to the translational invariance of the action, functional 
integrals in nonlinear field theory are not well defined when expanded about the 
classical solutions of the field equations. The translational symmetry results in zero 
eigenmodes of the second variation operator of the action which are a reflection of 
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the arbitrariness of locations in space of solitons or instantons (in Euclidean time). 
Obviously physical quantities are independent of such center-of-mass-like locations. 
A consequence of the existence of these normalizable "zero modes" which belong to 
the spectrum of linearised fluctuations in the soliton's background is a divergence 
one encounters when quantizing the theory in the background of these solitons or 
instantons. In essence the change of variables in field theory converts the zero modes 
into corresponding collective coordinates. However, calculation of the Jacobian of the 
transformation of-jvariables in the path integral formalism usually involves a Faddeev- 
Popov procedure™. 

Here we report an alternative method to carry out the transformation of variables 
and the evaluation of the related Jacobian in a systematic way with BRST invariance, 
and we apply it to quantum tunneling processes. The crucial point of the method 
is the exploitation of a shift invariance of the fluctuation action. As a test case 
the quantum tunneling effect in 1 + dimensions is calculated-iaiuthe Sine-Gordon 
(SG) potential. In the usual collective coordinate method 1131111113 time-dependent 
collective coordinates are associated with the D (spatial) dimensional (static) solitons 
indicating the motion of the kink or center of mass and one has to deal with phase- 
space functional integrals in 1 + D dimensions. Unlike this procedure our BRST 
invariant treatment does not evoke the time-dependence in one higher dimension 
and works on configuration-space path-integrals. Tunneling in the case of the Sine- 
Gordon potential is itself an interesting subject in view of similar features of the 0(3) 
nonlinear sigma model in the discussion towards understanding instanton induced 
baryon-number violating processes OJEJ. 



2. Energy band structure for the SG potential and the transition ampli- 
tude for quantum tunneling 

The 1 + dimensional Lagrangian we consider is 
with the SG potential 

y(0) = l[l + cos(#)] (2) 

where g > denotes a dimensionless coupling constant. The classical solution which 
extremizes the action is seen to satisfy the equation of motion which after one inte- 
gration with integration constant zero is 

i(f)Wno ,3) 

with euclidean time r = it. Mass m = 1 and natural units c = % = 1 are used 
throughout. The relevant classical solution which interpolates between neighbourung 
vacua, say, + = ^ and 0_ = — ~, is 

2 

<p c = - sin _1 [tanh(r + a)} (4) 
9 

where the integration constant a is interpreted as the position of the instanton. 
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In the following we consider the case of very high potential barriers with corre- 
spondingly small tunneling contributions to the eigenvalues. Thus we suppose |0)+, 
|0)_ are degenerate eigenstates in neighbouring wells with the same energy eigenvalue 

So such that H°\0)± = £o|0)± where H° is the Hamiltonian of the harmonic oscillator 
as the zeroth order approximation of the system. The degeneracy will be removed by 
the small tunneling effect which leads to the level splitting. The eigenstates of the 
Hamiltonian H then become 



|0)o 



1 



V2 



(|o>+-|o>. 



|0>, 



V2 



(|0>+ + |0>- 



(5) 



with eigenvalues S + AS and S 
oscillator level. It is obvious that 



AS respectively. AS denotes the shift of one 



.(0\H-H°\0)_ = AS. 



(6) 



These shifts embrace the energy bands which result from the translational invariance 
of V. We now calculate this energy resulting from tunneling with the instanton 
method. The amplitude for a transition between neighbouring vacua due to instanton 
tunneling can be written 



where 



(n — l)n 



S 



-T) = K((j)f, T; <j> h -T) 



V{<P}e 



-s 



(7) 




2 \dT 



+ V{<j>) 



eh 



is the Euclidean action. In the large time limit T — > oo, we have 

K = E^/l n >H e " 2flT l n/ ><nU> 

n,n' 

~ ^ o (0/)^o(</> l )e _2Tfo sinh(2A^T). 



(9) 



3. Translation invariance of action and zero modes 

As mentioned in the introduction the degeneracy of an action which possesses a 
translation symmetry leads to ill defined functional integrals in perturbation expan- 
sion about the classical configuration C since the symmetry results in zero modes of 
the second variation operator of the action. We consider this now in more detail. We 
expand 0(r) about the classical trajectory c (r) and so set 

0(r) = c (r) + X (r) (10) 

with the boundary conditions x(T) = x(~ T) = for the fluctuation. Substituting 
4>{t) of eq. flip]) into eq. (0) and retaining only terms up to the second order in \ f° r 
the one-loop approximation we obtain 

K = e~ Sc I (11) 
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where the classical action is evaluated along the trajectory 4> c so that 

2 



The fluctuation integral I is seen to be 




dr 



X(T)=0 
X(-T)=0 



V{x}e 



-AS 



with the fluctuation action 



where 



M 



AS 
2dr 2 ' 



T 



T 



1 



xMxdr 



cosh (t — a) 



(12) 



(13) 



(14) 



(15) 



is the self-adjoint operator of the second variation about the classical trajectory. 
Expanding x( T ) i n terms of normalized eigenfunctions of M we set 



X(r) 



(16) 



where 

M^ m = E m ^ m . (17) 

Changing the integration variables of flT3| ) to {c m }, the functional integral I can be 
formally evaluated to be 



dx(r> 



dCrr 



n 



7T 



18) 



which is seen to be divergent in view of the vanishing eigenvalue o 



is linear the Jacobian 



Ml) 



dc m 



the zero mode, 
is constant and 
y resorts to the 



Eq = 0. Since the transformation 

has therefore been factored out. TQ-,|t;ure the problem one normal 
so-called Faddeev-Popov procedureQo in order to transform the integral over the 
zero mode Cq into the continuous integration of a collective coordinate, which is the 
instanton position a in our case. 



4. BRST invariance and "gauge fixing" 

In the usual collective coordinate method, the essential ingredient is a change of 
variables which is such thaL-fiizery collective coordinate which is time-dependent is 
associated with a zero modeOll3. In our (1 + 0) dimensioiiaLcase the procedure is not 
appropriate. We therefore adopt an alternative methodli^Eil which employs a BRST 
invariance in dealing with the transformation of variables. If we identify euclidean 
time r with a spatial coordinate, the equivalence of the present method with that of 
a collective coordinate method with BRST invariance is similar to that demonstrated 
in a previous papenlB. 
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After expansion about the classical trajectory <; 
tains a shift symmetry expressed by the invariance 

AS( X ') = AS( X ) 



the fluctuation action still re- 



(19) 



where 



X = X + 



Oct 



x + 



dr 



(20) 



In other words the action is invariant under a kind of "gauge transformation" of 
the fluctuation variable x- This is an important observation. The key point of the 
BRST procedure is to enlarge the number of degrees of freedom and invent a nilpotent 
symmetry which mimics the structure of the gauge symmetry. Then one can achieve 
the effects of gauge fixing without breaking the BRST invariance of the system. To 
achieve this we first enlarge the configuration space by considering the parameter 
a as a variable. We replace the transformation fl2"U|) by the introduction of new 
anticommuting variables c and c and a Nakanishi-Lautrup auxiliary variable b (the 
latter in such a way that it implements the gauge fixing condition as its equation of 
motion, i. e. b = J X^dr in analogy with the implementation of the Lorentz gauge 
in QED) together with conjugate momenta P a , IT C , U 5 such that 



5x 

cki 
5c 
5b 
5c 



-c 



da 



—c 
2ixb 

0. 



(21) 



The variables a, c, c, x an d their conjugate momenta are assumed to satisfy the 
canonical Poisson relations 



{a, P a } 

{c,n c -} + 
{c,n c } + 
{x,k x } 



1 
1 
1 
1. 



(22) 



Then one finds that the BRST transformations can be generated by the following 
BRST charge 

d(p c 



-cP a - C7T 

da 



+ U- C b 



(23) 



which is nilpotent. The following BRST invariant term may now be added to the 
fluctuation Lagrangian to break the shift symmetry: 



27Tb 



da 
da 



■X 



dr — Tib 2 



Xdr 



cc 



d 2 4>c 
da 2 



X + 



'dtpc 
da 



dr — Tib 



(24) 
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Eliminating the Nakanishi-Lautrup auxiliary variable by using its equation of motion, 
we obtain the final expression of the fluctuation functional integral 



• cxp 



V{ X }V{c}V{c}V{a} 
xMx dr-TT 



da 



Xdr 



cc 



2 



da 2 



da 



dr } (25) 



(In the following we neglect the term linear in the fluctuation x since it is of higher 
order than the soliton mass M = J dr O^jjf) )• Expanding the fluctuation variable x 
in terms of the eigenmodes of M so that 

dx 



X 



dc r , 



T>{c m } 



(26) 



we obtain (ignoring the higher order contribution proportional to x in the coefficient 
of cc) 



dx 



dCn 



f V{c m }V{c}V{c}V{a} exp J - ]T c 2 m E m - nc 2 M 



ccM 



(27) 



where we used £0 = and the fact that \l/o is a normalized eigenfunction, i. e. (from 
(|) and O) 

*o = 4^, VM*„ = ^ (28) 



fM dr da 
where M = S c = is the classical action or instanton mass (as it is also called). 
Integrating out all variables (and recalling that T>{a} da) the final result is 

dx 



2T 



dc r , 



n 

m^O 



7T 



M 



= 2TI Vm 



where 



dx 



dc r , 



n(£ 



(29) 



(30) 



In / the factor comes from the T>{cq} integration and M from integrating out 

T>{c}V{c}. We see therefore that the functional integration can be done without re- 
sorting to the Faddeev-Popov method of inserting a delta function and interchanging 
integration and limiting procedures. The BRST procedure converts the ill defined 
integral over the zero mode into a Gaussian integral and leads to the integration over 
the instanton position, da, which gives rise to 2T. 



5. The one— instanton transition amplitude and the contributions from one 
instanton and an infinite number of instanton— antiinstanton pairs 

The following calculations are similar to those of the level splitting for the double- 
well potential^ and of the decay rate for the inverted double-well potentiala. The shift 
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method which is effectively defined by the transformation 

X(r) = £(r) 



(31) 



where £ = ^P- converts the functional integral into Gaussian form along with bound- 
ary constraints which are then implemented with a delta function-trick analogout to 
the Faddeev-Popov method which then leads to a finite result EJ. Taking care of 
the boundary condition constraint, the functional integral for the field fluctuations is 
given by 



Considering the large time limit one finds 



(32) 



T^oo 




(33) 



Comparing with eq. ([18]) it is seen that 

dx 



dcr, 



n 



/ 71 




(34) 



The eigenvalue Eq vanishes only asymptotically when time T tends to infinity; its fi- 
nite value corresponds to a finite-time interval. Eq can be evaluated with a so-called 
boundary perturbation methodOa which gives rise to a formula for the "unrenormal- 
ized" eigenvalue-zero eigenfunction £ 



dr ' 

d£(-T) 
dr 



En 



i 2 dr. 



(35) 



Again we take the large time limit; the quasi zero eigenvalue due to the size effect is 
obtained from eq. (|35| ) as 



Eq = 4e" 



-2T 



(36) 



Replacing J in eq. ( p9|) by eq. ( |34"D with the value Eq of eq. (0), the fluctuation part 
of the transition amplitude of the one-instanton sector is seen to be 



j(D 



2TVM— e 

7T 



(37) 



The contribution stemming from one instanton together with an instanton-paptiinstanton 
pair can be calculated with the help of the group property for propagators 113. The 



result is 



j(3) 



(2T) ; 
3! 



My 



V2 



TV 



(38) 
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where the determinant 



A 



7T 



1 d 2 S c 



(39) 



is determined from the end point integration of the group property of the propagator 
and is evaluated with a formula given in the literature^ by carefully taking the large 
time limit. The contribution from one instanton plus n pairs is a straightforward 
extension of eq. (|3~8D and is found to be 



j(2n+l) 



(2T) 



2n+l 



(2n + 1)! 



M) 



2n+l 




2ra+l 



A 2n e~ T . 



(40) 



The final result of the propagator is 



sinh 



i 

)5 \ 2 



2T 



2tt 



7 

e ^ 



(41) 



Comparing with eq. (|9|) we obtain the level shift, which is half of the width of the 
lowest energy band, i. e. 




A£ = J=-=[-\ (42) 



which is in exact agreement with the result in the literature^. The periodic potential 
has also been dealt with in the literature^ but the width of the energy band has not 
been given explicitly there. 
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